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Abstract —This paper experimentally and theoretically inves¬ 
tigates the fundamental bounds on radio localization precision 
of far-field Received Signal Strength (RSS) measurements. RSS 
measurements are proportional to power-flow measurements 
time-averaged over periods long compared to the coherence time 
of the radiation. Our experiments are performed in a novel 
localization setup using 2.4GHz quasi-monochromatic radiation, 
which corresponds to a mean wavelength of 12.5cm. We exper¬ 
imentally and theoretically show that RSS measurements are 
cross-correlated over a minimum distance that approaches the 
diffraction limit, which equals half the mean wavelength of the 
radiation. Our experiments show that measuring RSS beyond a 
sampling density of one sample per half the mean wavelength 
does not increase localization precision, as the Root-Mean- 
Squared-Error (RMSE) converges asymptotically to roughly 
half the mean wavelength. This adds to the evidence that the 
diffraction limit determines (1) the lower bound on localization 
precision and (2) the sampling density that provides optimal 
localization precision. We experimentally validate the theoretical 
relations between Fisher information, Cramer-Rao Lower Bound 
(CRLB) and uncertainty, where uncertainty is lower bounded by 
diffraction as derived from coherence and speckle theory. When 
we reconcile Fisher Information with diffraction, the CRLB 
matches the experimental results with an accuracy of 97-98%. 

Index Terms —Radio localization, Cramer-Rao Bounds, Fisher 
Information, Bienayme’s Theorem, Sampling theorem, Speckles, 
Uncertainty Principle. 


I. Introduction 

R ADIO localization involves the process of obtaining 
physical locations using radio signals. Radio signals are 
exchanged between radios with known and unknown positions. 
Radios at known positions are called reference radios. Radios 
at unknown positions are called blind radios. Localization of 
blind radios reduces to fitting these measured radio signals to 
appropriate propagation models. Propagation models express 
the distance between two radios as a function of the measured 
radio signals. These measured radio signals are often mod¬ 
eled as deterministic radio signals with noise using a large 
variety of empirical statistical models (5j. Radio localization 
usually involves non-linear numerical optimization techniques 
that fit parameters in the propagation model given the joint 
probability distribution of the ensemble of measured radio 
signals. Localization precision is usually expressed as the Root 
Mean Squared Error (RMSE). In the field of wireless sensor 
networks, the estimation bounds on localization precision are 
often calculated by the Cramer-Rao Lower Bound (CRLB) 
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from empirical signal models with independent noise GD- 

(IT Go). 

In general, localization precision depends on whether the 
measured radio signals contain phase information. Phase in¬ 
formation can only be retrieved from measurements that are 
instantaneous on a time scale that is short relative to the oscil¬ 
lation period of the signals. The smallest measurable position 
difference depends on how well phase can be resolved. This 
is usually limited by the speed and noise of the electronics 
of the system. Less complex and less expensive localization 
systems are based on measurements of time-averaged power 
flows or Received Signal Strengths (RSS). Time-averaging is 
usually performed over timespans that are large compared to 
the coherence time of the radios, so that phase information 
is lost. RSS localization is an example of such less-complex 
systems. Hence, when determining the bounds on localization 
precision, it makes sense to distinguish between the time scales 
of the signal measurements, i.e. between instantaneous and 
time-averaged signal and noise processing. 

Speckle theory describes the phenomenon of power- 
flow fluctuations due to random roughness from emitting 
surfaces. This theory shows that independent sources generate 
power-flow fluctuations in the far field that are always cross- 
correlated over a so-called spatial coherence region. The linear 
dimension of this region equals the correlation length. This 
correlation length has a lower bound in far-field radiation. 
In all practical cases of interest, the lower bound on the 
correlation length of this far-field coherence region is of the 
order of half the mean wavelength of the radiation. This lower 
bound is called the diffraction limit or Rayleigh criterion 0’ 
( 22 ). 

It is not obvious that this lower bound on correlation length 
holds in our wireless sensor network setup with its charac¬ 
teristic small cylindrical antennas. Therefore, we derive this 
lower bound on the correlation length with the corresponding 
cross-correlation function from the Maxwell equations using 
the IEEE formalism described by Our derivation of 

this lower bound on the correlation length and corresponding 
cross-correlation function leads to the well-known Van Cittert- 
Zernike theorem, known in other areas of signal processing 
like radar (29) , sonar (27) and optics (28) . Our novel experi¬ 
mental setup validates that power-flow fluctuations are cross- 
correlated over a correlation length of half a wavelength by 
increasing the density of power-flow measurements from one 
to 25 power-flow measurements per wavelength. 

In the field of wireless sensor networks, the estimation 
bounds on localization precision are usually determined by 
empirical signal models with independent noise over [space, 
time] fig, Q7), (20). Hence, the correlation length of the 
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radiation is assumed to be infinitely small. This representation 
renders practical relevance as long as the inverse of the 
sampling rate is large compared to the correlation length, so 
that correlations between measurements are negligibly small. 

This paper experimentally and theoretically investigates 
radio localization in the far field when sufficient measurements 
are available to reveal the bounds on localization precision. We 
use Bienayme’s theorem |9j §2.14] to show that the ensemble 
of correlated power-flow measurements has an upper bound 
on the finite number of independent measurements over a 
finite measuring range. We show that this finite number of 
independent measurements depends on the correlation length. 
When we account for the finite number of independent mea¬ 
surements in the Fisher Information, the CRLB on localization 
precision deviates 2-3% from our experimental results. We 
show that this approach provides practically identical results 
as the CRLB for signals with correlated Gaussian noise j6j, 
§3.9], given the lower bound on correlation length. Hence, 
the lower bound on correlation length determines the upper 
bound on localization precision. There are a few papers that 
assume spatially cross-correlated Gaussian noise on power- 
flow measurements m- These papers consider cross¬ 
correlations caused by shadowing that extend over many wave¬ 
lengths. p~8] states that their cross-correlation functions do not 
satisfy the propagation equations. Such cross-correlations have 
no relation to wavelength of the carrier waves as the diffraction 
limit is not embedded in the propagation model like we derive. 

Correlation, coherence and speckle properties of far-field 
radiation are governed by the spreads of Fourier pairs of 
wave variables that show up as variables in the classical 
propagation equations of electromagnetic waves. The products 
of these spreads of Fourier pairs express uncertainty relations 
in quantum mechanics, which are formulated as bandwidth 
relations in classical mechanics 03 and signal processing ||8j. 
0 mathematically establishes a relationship between Fisher 
information, the CRLB and uncertainty. As uncertainty is 
lower bounded by diffraction, Fisher information is upper 
bounded and CRLB is lower bounded. Our experimental 
work reveals quantitative evidence for the validation of this 
theoretical work. 

This paper is organized as follows. Section [TI] describes the 
experimental setup. Section III derives the propagation and 
noise models from first principles for our experimental setup 
for each individual member of the ensemble of measurements. 
Section [TV] describes the signal model. Maximum Likelihood 
Estimator (MLE) and CRLB analysis for the ensemble of 
measurements using the results of Section III Section [V] 
presents the experimental results in terms of spatial corre¬ 
lations and upper bounds on localization precision. Finally, 
Section [VT] provides a discussion and Section VII 
the conclusions. 


summarizes 


II. Experimental Setup 

Fi g. IT] and |2| show the two-dimensional experimental setup. 
Fig. Ill shows a square of 3 x 3m 2 , which represents the 
localization surface. We distinguish between two types of 
radios in our measurement setup: one reference radio and 
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Fig. 1. Measurement setup: reference radio is successively placed at known 
positions (crosses) on the circumference of the localization setup (3 x 3m 2 ) 
to measure RSS to the blind radio. The blind radio is located at the origin. 
For illustrative purposes, this figure only shows 60 of the 2400 measurement 
positions. 



Fig. 2. Photo of measurement setup: the blind radio is located underneath 
the aluminum table; the reference radio is located on the straight line in the 
picture, equidistantly positioned 0.5cm apart. The inset on the right shows a 
close-up picture of the radio. 


one blind radio. The reference radio is successively placed at 
known positions and is used for estimating the position of the 
blind radio. The crosses represent the 2400 manual positions of 
the reference radio (x-[. y-\ ... 22400 j 2/2400) and are uniformly 
distributed along the circumference of the square (one position 
every half centimeter). 

Rather than placing a two-dimensional array of reference 
radios inside the 3 x 3m 2 square, it may suffice to place 
a significantly smaller number on the circumference of the 
rectangle and get similar localization performance. Measuring 
field amplitudes on a circumference rather than measuring 
across a surface is theoretically expressed by Green’s theorem 
as is shown in Section III Whether sampling power flows on 
circumferences instead of sampling across two-dimensional 
surfaces suffices has yet to be verified by experiment. This 
paper aims to show the practical feasibility of this novel 
technique, which was first proposed by |23| §5.6]. 

The red circle represents the position of the blind radio. We 
place the blind radio at an unsymmetrical position, namely 
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(0, 0). We only use one blind radio and one reference radio to 
minimize influence of hardware differences between radios. 
The blind and reference radios are both main powered to 
minimize voltage fluctuations. The blind radio has a power 
amplifier and broadcasts messages with maximum power 
allowed by ETSI to maximize SNR. Both blind and reference 
radios have an external dipole antenna. The antennas have the 
same vertical orientation and are in Line-Of-Sight (LOS) for 
best reception. This implies that the polarization is vertically 
oriented perpendicular to the two-dimensional localization 
plane. The length of the antenna is half a wavelength. Its 
diameter is roughly one twentieth of a wavelength. We keep 
the relative direction of the printed circuit boards on which the 
antennas are mounted constant by realigning them every 25cm. 
In order to minimize interference from ground reflections, we 
place the radios directly on the ground, so that their antennas 
are within one wavelength height. The ground floor consists 
of a reinforced concrete floor covered by industrial vinyl. We 
minimize interference from ceiling reflections by placing a 
50 x 50cm 2 aluminum plate one centimeter above the blind 
radio antenna. The ground and the aluminum plate minimize 
the influence from signals in the z-direction, so that we only 
have to consider signals in two dimensions. All reference radio 
positions are in the far field of the blind radio. A photo of our 
setup is shown in Fig. [2] 

At each of the 2400 reference radio positions, the reference 
and blind radio perform a measurement round. A measure¬ 
ment round consists of 500 x 50 repetitive multiplexed RSS 
measurements to investigate and quantify the measurement 
noise and apply CRLB analysis to this measurement noise. 
Each measurement round consists of 50 measurement sets that 
consist of 500 RSS measurements on an unmodulated carrier 
transmitted by the blind radio. Between each measurement 
set of 500 RSS measurements, the reference and blind radios 
automatically turn on and off (recalibrate radios). Although we 
did not expect to find any difference from these two different 
forms of multiplexing RSS signals, our experiments should 
verify this, which they did. 

Hence, a measurement round consists of 50 measurement 
sets, each set consisting of 500 repetitive multiplexed RSS 
measurements. Measurement rounds and measurement sets 
are represented by = -Pi,i,i • • • ^ 2400 , 50 , 500 - Index l 

identifies the position of the reference radio and thus the mea¬ 
surement round, index m identifies the 50 measurement sets, 
and index n identifies the 500 individual RSS measurements in 
each measurement set. Pi m n represents the measured power 
in dBm and is time-averaged over 125ps according to the radio 
chip specification ED- The coherence time is 25ps, which 
implies a coherence length of 7.5km. The averaging time of 
125ps is a factor of five larger than the specified coherence 
time of the carrier wave of a typical 802.15.4 radio Eil- 
Practically, this means that these power measurements lose 
all phase information \26 §1.5]. In theory, this means that we 


measure the time-averaged power flow or Poynting vector as 
the cross-sections of the antennas are the same. 

The blind radio transmits in IEEE 802.15.4 channel 15 in 
order to minimize interference with Wi-Fi channels 1 and 6. 
The reference radio performs power-flow measurements in the 


same channel and sends the raw data to a laptop over USB, 
which logs the data. We use Matlab to analyze the logged data. 
Between each measurement round, we change the position of 
the reference radio by 0.5cm and push a button to start a new 
measurement round. Note that 0.5cm is well within the A/2 
diffraction limit of half the mean wavelength. 

In summary, each measurement set takes one second; each 
measurement round takes roughly one minute (50 x 1 seconds). 
The experiment consists of 2400 measurement rounds and 
takes in total 40 hours (2400 minutes « 40 hours). In practice, 
we spend almost three weeks in throughput time to perform 
this complete data collection. 


III. Propagation and Noise Model 
This section formulates the propagation and noise model 
of radio localization systems operating in the far field. This 
model holds for each individual member of the ensemble 
of 2400 power-flow measurements over space described in 
Section [U] Our propagation model is based upon EM inverse 
source theory, where we follow the IEEE formalism given 
by 1241. This formalism is needed to derive the lower bound 
on correlation length of wireless networks using power-flow 
measurements on a time scale long compared to the coherence 
length of the hardware. 


A. Measurement Configuration 


The measurement configuration of practical interest is com¬ 
posed of the wave-carrying medium, i.e. the atmosphere, a 
number of source domains and one receiver domain. Our 
measuring chamber has a variety of contrasting obstacles like 
a ceiling, walls, a ground floor, pillars, and a small aluminum 
table hiding the ceiling from the transmitter. Most of these 
secondary sources are so far away from the receiver that they 
are neglected in our theoretical representation. We account for 
the following domains; 


the wave-carrying medium denoted by R 3+1 , with electric 
permittivity eo > 0, magnetic permeability /iq > 0 
and EM wave speed Co = (eo/To) -1 ^ 2 . It is locally 
reacting, spatially invariant, time invariant and lossless. 
To locate positions in space, the observer employs the 
ordered sequence of Cartesian coordinates {x, y. z \ £ 
l 3 , or a; £ R 3 , with respect to a given origin O , 
while distances are measured through the Euclidean norm 
\x\ = sjx 2 + y 2 + z 2 . The fourth dimension is time and 
is denoted by t. 

a Transmitter with bounded support D £ R 3 with 
D = 1) J IJ D K . Here, l) J denotes the spatial support car¬ 
rying electric currents with volume density .//., and 1) K 
denotes the spatial support carrying magnetic currents 
with volume density \Ki J \~.The transmitter transmits 
unmodulated EM currents at a temporal frequency of 
loq/2tt = 2.42GHz with a temporal frequency bandwidth 
of Aojo/ 27r = 40kHz 130 ]. The polarization of the EM 
currents is assumed to be perpendicular to the ground 
plane of localization space, which is assumed to be the 
{x, y} plane. The electromagnetic volume current densi¬ 
ties are usually represented by equivalent surface current 
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densities, which physically relate to the bounded charges 
and currents in terms of polarization and magnetization 



• Volume scatterer(s) or volume noise with bounded spatial 
support D ' e R 3 with D v = D V J U D v - K . Here, 
l) V J denotes the spatial support carrying electric currents 
with volume density &J}', and !)'' J< denotes the spatial 
support carrying magnetic currents with volume density 
S [K}j] . We assume this noise to be negligible in our 
setup as we identify it with the thermal noise in all 
contrasting media in the measuring chamber. 

• Surface scatterer(s) or surface noise with bounded spatial 
support dD s e R 3 with dD s = dD s ’ J U dD s ’ K . Here, 
dl) s J denotes the surface boundary carrying electric 
currents with surface density 6J^, and dD s,R denotes 
the surface boundary carrying magnetic currents with 
surface density S [if? ■] . The surface noise is assumed 
to be caused by surface roughness of the ground floor as 

described by G3- 

• a Receiver with bounded spatial support D n in which 
the electric field strength E r and the magnetic field 
strength [H PA \ are accessible to measurement through 
the received time-averaged power flow 

P P =,l j T ""' 2 [H p , q ]-E q dt. ( 1 ) 

lm J—Tm/2 

The receiver measures the power flow of an unmodulated 
carrier with an observation time of T m = 125ps, long 
compared to the coherence time of 2ops (3g. 


B. Propagation Model 

The mapping SOURCES => FIELD is unique if the 
physical condition of causality is invoked. For the received 
signal at the point, x, of observation, the relation between 
the EM vector and scalar potentials, jAfc, j, of the 

received power flow and the real EM currents on the surfaces 
of the bounded sources is given by 124 §6] 


jAfcJT'jj] | (as,t) = G(x,t) ( * )( * | J fe , [Kij] }(a;,t), 

( 2 ) 

where the spatial and temporal convolutions extend over the 
surfaces of the spatial-temporal supports of the pertaining 
bounded primary and induced sources located at x D in the 
measurement configuration. According to Section |III-A| the 
spatial supports are the surfaces of the transmitter and of the 
ground floor. In two dimensions, the surfaces become circum¬ 
ferences of their cylindrical cross sections. Equation ([2} links 
the extended sources to the well-known retarded potentials. 
In our measurement configuration, where the polarization of 
transmitted and received signals are directed perpendicular to 
the {x,y} plane, the array of Green’s functions that couples 
the sources to their respective potentials takes on the scalar 
form 


Ait\x\ 


In (|3j, 8(t — \x\ /co) is the (3+l)-space-time Dirac distribution 
operative at x = 0 and t = 0 



/•Oo 

r. / 

MY 

/ exp 

zw 0 

t- — 

J —oo 

L V 

Co ) J 



(4) 


The bandwidth, Awo, of temporal frequency spectrum of the 
hardware is usually mathematically represented by its inverse, 
the coherence time, r c , as described by [ [26] §10.7] 

(87t ln2) 1 / 2 


Aw q 


(5) 


where we assume a narrow Gaussian line shape relative to 
the carrier frequency as is the case for quasi-monochromatic 
radiation 

4=5 « 1. (6) 

w 0 


Time and position are connected by the constant speed of light, 
cq, as are angular temporal frequency and wavelength Aq 


27TCo 27T 

Wo fc 0 


(7) 


To simplify the notation, we represent this time dependence 
of the source currents in the Dirac distribution by limiting the 
integration interval in 0 to the frequency bandwidth, Awq, 


t -= exp 

co 


zw 0 


Co 


sine Awo t — 


co 


Awq 
27r 


( 8 ) 


For Awo —> 0, the slowly varying envelope of the sinc- 
function approaches unity, and we end up with a plane- 
wave representation leading to a spherical wave in <(3j. The 
sinc-function can be replaced by a normalized Gaussian or 
Lorentzian as is applied in [26 §10.7]. For our work, the 


only time dependence that is important is that our observation 
time is long compared to the coherence time so that all phase 
information is lost. 


C. Far-Field Approximation 

The convolution integral over the extended spatial supports 
in (|2) assumes its far-field approximation when the point of 
observation, x, is far away from any part of the extended 
spatial supports, the coordinates of which are denoted by x D , 
so that 

-^ + o(M), (9) 

where x denotes the unit vector in the direction of observation, 
x. In the far-field approximation, the potentials in 0 are 
replaced by the far-field electromagnetic field strengths [24] 
§6] with the given polarization along the z-axis 

E™(x, t) = G°°(x : t) ( * K * J z (x,t), (10) 

where the tilde, ~ denotes the sum of a macroscopic average 
of all primary and secondary fields including a small noise 
term due to the surface roughness of the ground floor 

Ef{x,t) = Ef {x, t ) + SE™(x, t). 


\x\ \x , or \x — x \ = \x\ 


G(x, t) 


(3) 


( 11 ) 
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In ( flOl ), the far-field Green’s function is given by substituting 
( [Tol l in 0 and ([ 8 ]) in <[4]). The transverse far-field mag¬ 


netic field strength, H£°(x,t) = H£° z (x,t) , is directed 

perpendicular to the z-axis and x and equals the far-field 
electric field strength multiplied by co£o- In the far-field 
approximation, the time-averaged power flow of 0 becomes 
inversely proportional to |cc | -2 with the help of 0 , so that 
without any noise, 0 reduces to 


P{x,Tm Tc) — aj 0 1 


*0 


X 


2 ' 


( 12 ) 


In ([12]), P denotes the absolute value of the time-averaged 
Poynting vector given by 0 . and P\x 0 \ denotes the reference 
power flow at reference far-field distance of, say, |aJo| = lm. 
Equation © can be read as the non-logarithmic gain model. 


D. The Inverse-Source Problem 

In the inverse-source problem formulation pertaining to our 
measurement configuration, information is extracted from an 
ensemble of measured, time-averaged power-flow signals. This 
information reveals the nature and location of the scattering 
volume and surface sources in ([2]). The mapping of the 

Observed Field Ef°(x,t) =>■ Generating Sources J z (x D ) 

(13) 

is known to be non-unique. A detailed analysis in this respect 
can be found in ©■ In radio localization, an algorithm is used 
that is expected to lead to results with a reasonable degree 
of confidence. Such an algorithm is based on the iterative 
minimization of the norm of the mismatch in the response 
between an assumed propagation model for each individual 
measurement and an ensemble of observed power-flow signals. 


E. The Influence of Noise 

The influence of noise can be accounted for by using 
an input power-flow signal perturbed by an additive noise 
signal as is usually done in the Log-Normal Shadowing Model 
(LNSM). For independent, uncorrelated noise, © reduces to 
the non-logarithmic gain model |5] 

P(x, T m » T C ) = P M (lO*/ 10 ) , 

for E [A ' 2 ] 1/2 /10 < 1. (14) 

In ( fi~4| , X denotes the Gaussian variable of the independent 
zero-mean noise with standard deviation E [A' 2 ] 1 ^ 2 . Under the 
condition that the perturbation is small, we obtain 

10 x/1 ° ^ (1 + 2.83A/10), for E [X 2 ] 1/2 /10 < 1. (15) 

The linearization of © is established by recursive expansion. 
With this linearization of the perturbed, time-averaged power 
flow with added noise, the nature of the correlations can be 
derived from (JTj). ([0, ([0 and 0 and then be combined with 
( fl4| > and ( p~5j ). We apply the rationale used by Van Cittert and 
Zernike. 

According to this rationale, each surface element on the 
contrasting surfaces acts as a spatially incoherent point source 


represented by the Green’s function of 0. The contrasting 
surfaces observed by the receiver in the far-field in our local¬ 
ization setup are the surfaces of the primary and secondary 
sources. The primary source is the blind transmitter, and the 
main secondary source in our measuring chamber is the ground 
floor. The wave vector of the primary source is fixed in size, 
i.e. k()X, as the primary source is not extended and can be 
considered as a point source. The wave vectors k 0 x of the 
plane wave amplitudes originating from secondary stochastic 
point sources induced under grazing incidence all over the 
ground floor are the same in size. All these scattered plane 
waves that are directed towards the vertical antenna of the 
receiver are again directed parallel to the ground floor. As 
the vertical receiving antenna has a circular cross-section, the 
projections of all these wave vectors on the planes tangential 
to this cross section give rise to a spatial-frequency bandwidth 
of o k < k 0 . 

The expectation value, E, of the superposition of these plane 
waves in the far field at a slightly displaced position, x + Sx, 
all of the same amplitude but with random phases follows 
from 0, ([3]) and 0 


E[5E™ {x + 8x)\ = E 


(x-\-8x) 

G°°(x + 6x) * 6J‘(x) 


= E [exp [iko x ■ Sx] 5Ef ) (x)] 

I Sx I 

= E [exp [ikox ■ Arc]] E [SEf (a:)], with p—p <C 1, (16) 

|a;| 


where we have taken out the time dependence to simplify 
the notation. In ©, the third equal sign results from the 
spatial invariance of the propagating medium and from the fact 
that the variance in vertical surface roughness is statistically 
independent of horizontal correlations in the far-field. The 
ensemble average of all individual measurements over space 
of exp [ikgx ■ 6x] can be replaced by an average over the 
Fourier conjugated space because of the assumed wide-sense 
stationary (WSS) random noise process |28] §3.4], In one 
dimension, the expectation value of exp [ik^x ■ Sx] equals the 
sinc-function 


1 f k ° 

E [exp [ik 0 x ■ Acc]] = —— / exp [ikx ■ Arc] dk 

= sinc(k 0 Sx), (17) 


and in two dimensions the jinc-function jinc(fcoAcc), with 
5x = \5x\. The first-order spatial correlation function of the 
perturbed electric field strengths follows from ( fl 6 | ) by mul¬ 
tiplying its deterministic form for Ef > {x + Sx ) by E*’°°(x) 
and by neglecting the second-order terms of the perturbed field 
strength. We then take the expectation value 


E 

Ef°(x + Sx)E* 

00 (x) 

E 

'\E?(x)f 



= 1 + 


E 

5\E?(x)\ 2 


E 

\E?(x)\ 2 



sinc(fc 0 <5:r), 


(18) 

where the superscript star, *, denotes complex conjugate and 
the ratio in the right member denotes the normalized variance 
of the perturbed field. 

The sine- and jinc-functions have their first zeros at Sx = 
Aq/ 2 and Sx = O. 6 IA 0 , which are called the diffraction limit 
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or Rayleigh criterion in one and two dimensions. They can 
be considered as the far-held correlation functions for each 
individual measurement with a spatial correlation length that 
equals the diffraction limit. Equation ( fl8| > is equivalent with 
the Van Cittert-Zernike theorem (22). This theorem links far- 
held spatial correlations of noise to diffraction. 

The diffraction limit acts as a fundamental lower bound 
on the spatial correlation or coherence region for time scales 
Tm t c , and hence, as an upper bound on the localization 
precision as we will see in Section IV Equations ( fl6l > - ( fT8j ) 
generally allow for extracting the phases from the EM held 
measurements as long as the timescales are short relative to 
an oscillation period T m <C 2tt/ u>q. Under those time scales, 
the Van Cittert-Zernike theorem is equivalent with the Wiener- 
Khintchine theorem. 

Using Q and ©> we extend the perturbed input signals 
of ( p~4l > and ( | 1 5[ i from independent to second-order spatially 
correlated noise to 


E [P(x)P{x + (Sa;)] 
E [P 2 {x)] 


= 1 + 


E 

6\E?(x)\ 4 ' 

E 

'\E?(x)\* 


\k 0 Sx) 


E Tx 2 1 sine 2 ( fco <5 x ) n r 

S* 10-ro-, for E [X 2 ] /TO < 1, T m » 


T r , 


(19) 

The normalized scattering cross-sections in © and © 
resulting from surface roughness equal the measurement vari¬ 
ance E [X 2 ] at each individual reference-radio position. As 
we shall see in the next subsection, this measurement variance 
has a fundamental lower bound within the time scale of our 
measurements ( T m >> t c ). Although the correlation function 
derived only holds for correlations in the close neighborhood 
of each point of observation, © can be extended to larger 
distances when the path losses are accounted for as is usually 
done in LNSM. 

The correlation function adopted in the literature is an 
exponentially decaying function |4), [18). Reference )T8) notes 
that the exponentially decaying correlation function does not 
satisfy the propagation equations of EM theory. Our correla¬ 
tion functions hold for each individual measurement position 
of the reference radio. The extension to an ensemble of 
measurement positions with a derivation of the measurement 


variance of this ensemble is given in Section IV Section [TV] 
investigates two signal models for this extension and compares 
the computed measurement variances for the two correlation 
functions of power flows as a function of the sampling rate. 


F. Uncertainty and Fundamental Bounds 

Equation © contain two products of the Fourier pairs 
[angular frequency w, time t] and [spatial frequency k, position 
a;]. These two sets of conjugate wave variables naturally lead 
to an uncertainty in physics called the Heisenberg Uncertainty 
principle. Hence, uncertainty is a basic property of Fourier 
transform pairs as described by © §2.4.3] 

overt > 0.5, (20) 

and 

( 21 ) 


Equations ( |20| ) and ( |2T| are the well-known uncertainty re¬ 
lations for classical wave variables of the Fourier pairs [en¬ 
ergy, time] and [momentum, position]. In these Fourier pairs, 
energy and momentum are proportional to angular frequency 
and wave number, the constant of proportionality being the 
reduced Planck constant. Although these uncertainty relations 
originally stem from quantum mechanics, they fully hold in 
classical mechanics as bandwidth relations. In information 
theory, they are usually given in the time domain and are 
called bandwidth measurement relations as described by |8). 
In Fourier optics, they are usually given in the space domain 
as described by ©. 

Heuristically, one would expect the lower bounds on uncer¬ 
tainty in time, a t , and position, a x , to correspond to half an 
oscillation cycle, as half a cycle is the lower bound on the 
period of energy exchange between free-space radiation and a 
receiving or transmitting antenna giving a stable time-averaged 
Poynting vector. When one defines localization performance 
or precision as the inverse of the RMSE, one would expect 
the localization precision of the ensemble not to become 
infinitely large by continuously adding measurements. When 
the density of measurements crosses the spatial domain of 
spatial coherence as computed from ©, the power-flow 
measurements become mutually dependent. The link between 
uncertainties and Fisher Information was first derived by |2). 
As we show in Section [V] our experiments reveal convincing 
evidence for this link. 


IV. Signal Model, Maximum Likelihood Estimator 
and Cramer-Rao Lower Bound 

The first five subsections describe the signal model, MLE 
and CRLB usually applied in the field of RSS-based ra¬ 
dio localization for independent noise ©’ EH©- © 
and cross-correlated noise 0, ©■ Section |IV-F| reconciles 
diffraction and the Nyquist sampling theorem with the CRLB 
using cross-correlated noise and Bienayme’s theorem. In ad¬ 
dition, it introduces the notion of the maximum number of 
independent RSS measurements and relates this to Fisher 


Information. In Section IV-G simulations verify that the 
estimator is unbiased and efficient in the setup as described 
in Section m We extend the notations introduced in Section 
|H| and III that describe our measurement configuration with 
the bold-faced signal and estimator vectors usually employed 
in signal processing. The experimental setup consists of a 
reference radio that measures power at 2400 positions xi = 
( xi,yi ) G .. .a;24oo,?/24oo} of an unmodulated carrier 

transmitted by the blind radio. At each position, the reference 
radio performs 50 x 500 = 25,000 repetitive multiplexed 


— -Pi, 1,1 • • ■ P2400, 


50,500 


. These 


power measurements Pi, 
measurements are used to estimate the blind radio position, 
which is located at the origin x D = ( x,y) = (0,0). 


A. Empirical Propagation Model 

We adopt the empirical LNSM for modeling the power 
over distance decay of our RSS measurements. As the cross- 
sections of the blind transmitter and reference receiver are 


over > 0.5. 
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given and equal, the power as well as the power-flow mea¬ 
surements are assumed to satisfy the empirical LNSM |5) 

Pi = P(ri) + Xi, ( 22 ) 

where 

P(n) = Pr 0 - 10 r/ log 10 ^ . (23) 

In ( [22] ) and ( |23[ ), i identifies the power-flow measurement. 
P(rj) denotes the ensemble mean of power-flow measure¬ 
ments at far-held distance 


n = \Xi - x 


D | 


(24) 


in dBm. P ro represents the power flow at reference distance 
ro in dBm. 77 represents the path-loss exponent. X t represents 
the noise of the model in dB due to fading effects. X,; follows 
a zero-mean Gaussian distribution with variance cr dB and is 
invariant with distance. Equations ( |22| ) and ( |23| are equivalent 
with the 10 log 10 of where the path-loss exponent equals 
two. 

One usually assumes spatially independent Gaussian noise 
(H), Gz)’ GD’ HD- in GD’ the cross-correlation between 
power-flow measurements is independent of wavelength and is 
modeled with an exponentially decaying function over space 
by 

p(5x ) = exp ^ 


L X J 


(25) 


where x denotes the correlation length. Section 1II-E shows 


that the lower bound on correlation length depends on the 
wavelength and equals half the mean wavelength, X — A 0 /2. 
In Section |III-E| ( fl9| shows that the cross-correlation func¬ 
tion of power-flow measurements takes on the form of the 
diffraction pattern 0 - 0 


p(Sx) = sinc 2 (fcofe). 


(26) 


In Section IV-E we show the implications of using the cross¬ 


correlation function of ( |25) > and the cross-correlation function 
of ( |26| ) that satisfies the propagation model derived from first 
principles. 


Here, P(r r ) denotes power flow at far-field distance r, as 
expressed by ([24]). In ( |27[ i, G denotes the set of unknown 
parameters, where 6 C [x, y. P ro . 77 ]. C denotes the covariance 
matrix. Usually power-flow measurements are assumed to 
be independent so that C = where I is the identity 

matrix p2] , (17), (19), |2T]|. In case of correlated noise (4), 
0 , the elements in the covariance matrix are defined by 
Cij = Pip<J 2 dB , where 


Pip = p(n,j), (31) 

where 

r i,j = \xi-Xj\ (32) 

denotes the correlation distance between reference radio posi¬ 
tions i and j. Equation ( |3 1 [ 1 denotes the correlation coefficients 
as expressed by ([25]) and ([26]). 


C. Calibrate Log-Normal Shadowing Model 

Usually, the three parameters [P ro , 77 , er^s] of the LNSM 
are calibrated for a given localization setup jT2), Q7), (21). 
The blind radio position is assumed to be known when the 
LNSM is calibrated. The LNSM assumes that OdB is equal for 
each RSS measurement, so that P ro and 77 can be estimated 
independent of o^b 

[P r cal ,r 7 cal ] =arg min V(6). (33) 

9=[ P r 0 M 

In ( |22j i, a d j\ is defined as the standard deviation between 
the measurements and the fitted LNSM using and // Kal . 
Equation ( [33] ) gives for all our power-flow measurements the 
best fit when the LNSM parameters are calibrated at 

P r cal = —16.7dBm, ? 7 cal = 3.36, = 1.68dB. (34) 

We use these calibrated LNSM parameter values to calculate 
the bias and efficiency of our estimator from simulations. In 
addition, we use these calibrated LNSM parameter values to 
calculate the cross-correlations in the far-field and CRLB. 


B. General Framework Log-Normal Shadowing Model 

The non-linear least squares problem assuming the LNSM 
with correlated noise is denoted by 

argminU(S) = argrnin (P — P(r)) T C -1 (P — P(r)), 

0 0 

(27) 

where P denotes the vector of power-flow measurements 

P = [P 1 ,...,P n ], (28) 

where P(r) denotes the vector of power-flow measurements 
calculated by the LNSM as expressed by ( |23j ) 

P= [P(n),...,P(r n )], (29) 

where r denotes the vector of distances between the reference 
radio and the blind radio 

r = [r 1 ,...,r n ]. (30) 


D. Maximum Likelihood Estimator 

We use the MLE as proposed by pT) , with the phys¬ 
ical condition of causality invoked on 77 and with boundary 
conditions on (x,y), to estimate the blind radio position 


[x mle , t/ mle , P“ le , r; mle ] = arg min V(6) 

9=[x,y,P T0 ,r)\ 

subject to 

77 > 0 

- 1 < x < 2 

- 2 < y < 1 , 


(35) 


Our estimator proces ses the noise as independent, so that C = 
Icr dB - Section [fV-G 


shows that this estimator is unbiased and 
efficient in our simulation environments with uncorrelated and 
correlated noise. 












E. Cramer-Rao Lower Bound 


CRLB analysis provides a lower bound on the spreads 
of unbiased estimators of unknown deterministic parameters. 
This lower bound implies that the covariance of any unbiased 
estimator, COV (G), is bounded from below by the inverse of 
the Fisher Information Matrix (FIM) F as given by |6j| 

CO W(G)>F~ 1 (G). (36) 


In ( |36l >, G represents the set of unknown parameters, and 
G represents the estimator of these parameters. In case of 
multivariate Gaussian distributions, the elements of the FIM 
are given by @ §3.9] 


imu = 


'dP' 

T 

C - 1 

<9P' 

[d0 a \ 




r -i9C ,dC 

dd a dO b 


(37) 


In our case, the covariance matrix is not a function of unknown 
parameter set G so that the second term equals zero. The 
elements of the 4x4 FIM associated with the estimator 
defined by [35) , G = [x, y. P ro , 77 ], are given in the independent 
measurement case by 03 and in the correlated measurement 
case without nuisance parameters by tH- The elements of F 
consists of all permutations in set G. After some algebra, 
reduces to 


<9P 

dxi 

dP 

dyi 

dP 

dr/i 

dP 

dPZ, 


= -b 


{X - Xi ) 


= -b 


( y~Vi ) 


10 

ln( 10 ) 
= 1 , 


ln(rj) 


where 


We use 




(38) 


(39) 


to calculate the CRLB of the MLE expressed 
by ( [35]) using the calibrated LNSM parameter values given 
by ({34ji: \a dB = cr c d B,r] = r/ cal ]. The CRLB on RMSE for 
unbiased estimators is computed from 


RMSE = W / E((a; mIe - x) 2 + (y m,e - y ) 2 ) > ^/tifF-^G)). 

(40) 

Here tr(-) represents the trace of the matrix. 

Fig. [3] shows the lower bound on the RMSE calculated 
by [38|l and [40) as a function of the number of RSS 
measurements. When we assume independent RSS measure¬ 
ments, C = Icr db , the RMSE decreases to zero with an 
ever increasing number of independent RSS measurements. 
This is in accordance with the theoretical bound analysis 
presented by |l2), [17], [19], |20) . Hence, there is no bound 
on localization precision with increasing sampling rate. When 
we account for the spatial cross-correlations between RSS 
measurements, Cij = Pijcr dB , the bound on radio localization 
precision reveals itself when sufficient RSS measurements are 



Fig. 3. Theoretical RMSE calculated by CRLB as a function of the number 
of RSS measurements over space. The green curve shows the RMSE as a 
function of the number of independent RSS measurements over space. The 
blue curve shows the RMSE as a function of the number of RSS measurements 
with cross-correlations equal to the diffraction pattern as expressed by (26) . 
The blue triangles show the RMSE calculated by the corresponding CRLB 
as expressed by (38) . The covariance matrix becomes ill-conditioned when 
measurement density goes above a certain threshold (see text). The black 
curve shows the RMSE as a function of the number of RSS measurements 
with an exponentially decaying cross-correlation as expressed by (25) . The 
black triangles show the RMSE calculated by the corresponding CRLB as 
expressed by (38) . These results are obtained using the same setup as described 
in Section |TT] We use the calibrated LNSM parameters as given by El- 


available. We evaluate the CRLB with the cross-correlation 
function expressed by the diffraction pattern in [26) and the 
exponentially decaying cross-correlation function as expressed 
by ( |25| >. They basically converge to the same bound on 
localization precision. Apparently, the bound on localization 
precision depends on the bound on the correlation length rather 
than on the form of the cross-correlation functions. 

The determinant of the covariance matrix starts to decrease 
when correlations start to increase, until the Fisher Information 
and thus localization precision stabilizes on a certain value. 
When the cross-correlations are equal to the diffraction pattern 
as expressed by [26) , the covariance matrix becomes singular 
when the sampling rate goes above a certain threshold, which 
equals the inverse of the diffraction limit or the Nyquist 
sampling rate. Hence, the multivariate Gaussian distribution 
becomes degenerate and the CRLB cannot be computed with¬ 
out regularization. 

F. Diffraction, Sampling Theorem, Covariance and Fisher 
Information 

This section connects diffraction to the sampling theorem, 
covariance and Fisher Information using Bienayme’s theorem 
|9j §2.14], Bienayme’s theorem offers a statistical technique 
to estimate measurement variance of correlated signals. It 
states that when an estimator is a linear combination of n 
measurements X - t ... X n 

n 

0 = ^w i X i , (41) 

2 = 1 
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the covariance of this linear combination of measurements 
equals 


CO v(d) = J2 w * a *+Y,J2 WiWjPijcriCTj (42) 

i—1 i= 1 

and is equivalent to the measurement variance. In ( [42] ), Wi 
is a weighting factor, <Ji represents the standard deviation of 
measurement i, and pij represents the correlation coefficient 
between measurement Xj and X ;) . When all measurements 
have equal variance of = of = a 2 and equal weights Wi = 
Wj = 1/n, ( |42| ) reduces to 

COV(0) = (— + —— a 2 = -^—a 2 . (43) 

\n n ) riefr 

In (|43), p is the spatially averaged correlation of n[n — 1) 
measurement pairs and is given by 

^ n n 

P = 7-tv tH -:i (44) 

n[n — 1) ^ 

' *=i 


and n e n represents the number of measurements that effec¬ 
tively decreases estimator covariance and is given by 


n e ff = 


1 + (n— l)p 


(45) 


0 shows in the time domain that the cross-correlation of 
bandwidth limited fluctuations on signals with a periodic wave 
character is proportional to the Point-Spread-Function (PSF) 
of a measurement point. In the space domain, the cross¬ 
correlation of far-field power-flow measurements takes on the 
form expressed by the diffraction pattern in ( |26| . In addition, 
|3j shows that the radius, 5rij, of the first zero 


p(r itj = A 0 /2) = 0 (46) 


equals the inverse of the minimum (Nyquist) sampling rate. 
The sampling theorem determines the sampling rate that 
captures all information of a signal with finite bandwidth. 
For signals with a periodic wave character measured in the 
far field, noise is cross-correlated over the far-field spatial 
coherence region of the signal. This far-field coherence region 
assumes the far-field diffraction pattern and is equivalent to 
the PSF of power-flow measurements Q, 0 . 

The lower bound on the far-field region of spatially cor¬ 
related signals and noise is equal to the diffraction limit 
0 - © and to the inverse of the upper bound on spatial 
frequencies divided by 2n, i.e. (2k’o/27r) _1 . The finite spatial- 
frequency bandwidth of the far field relates diffraction to the 
sampling theorem 0 and Whittaker-Shannon interpolation 
formula or Cardinal series [14, §2.4]. Hence, this coher¬ 
ence region determines the maximum number of independent 
measurements over a finite measuring range. This maximum 
number of independent measurements is also known as the 
degrees of freedo m |l] , and is given by ?r e ff in ( |45j ). 

We decompose (|43| into two factors 


COV(0) 



(47) 


where ( ff J represents the ratio of the total number of 
measurements and effective number of measurements, and 


represents the covariance of the estimator assuming that 
all measurements are independent. The global spatial average 
of correlation coefficients, p, as given by ( |44[ > can be computed 
using the correlation function ( |26] > for any localization setup by 
assuming an ever increasing sampling rate. For our localization 
setup of Fig. |T] we compute p k, 0.0048 processing all 2400 
RSS measurements over space. Substituting this value in ( |45j ) 
gives n e g ~ 191. 

To account for the effects of spatial correlations on the 
covariance of non-linear unbiased estimators, we heuristically 
rewrite ( |36] > into the following postulate 

COV(0) > (— ) F~ l {G), (48) 


where F is equal to the Fisher Information assuming indepen¬ 
dent measurements. Our postulate rigorously holds when the 
variance and thus the Fisher Information of each wave variable 
at each measurement point is equal, assuming that the CRLB 
holds as defined in ( |36| . Our postulate of ( |48j ) is in agreement 
with the concept that the degrees of freedom of signals with 
a periodic wave character is independent of the estimator. 
Postulate ( |48] i is also in line with that CRLB efficiency is 
approximately maintained over nonlinear transformations if 
the data record is large enough [ 6 ] §3.6]. Hence, ( [48] ) is a good 
performance indicator when the correlation length is small 
relative to the linear dimensions of localization space. 

Fig. [3] shows the lower bound on the RMSE calculated 
by ([4~8p and ( [40] ) using the two cross-correlation functions 
expressed by ( [25] ) and ( |26] >. The solid curves represent the 
results of ( |48| >. Our approach expressed in ( |48j ) differs less 
than 1mm from the CRLB for signals with correlated Gaus¬ 
sian noise [ 6 ] §3.9] over the entire range shown by Fig. [3] 
The influence of correlated measurements on the CRLB is 
apparent for n > n e g as ( |48] i converges asymptotically to 
roughly half the mean wavelength for both cross-correlation 
functions. Hence, the correlation length and thus the finite 
number of independent RSS measurements determines the 
bound on localization precision. Note that the effective number 
of measurements converges to practically identical values for 
both cross-correlation functions, n e g ~ 191. The sine-squared 
cross-correlation function expressed by ( [26] ) is bandwidth 
limited, so that the covariance matrix becomes degenerative at 
roughly the Nyquist sampling rate. The exponentially decaying 
cross-correlation function expressed by ( [25] ) is not bandwidth 
limited, so that the covariance matrix stays full rank. In 
this case, the Fisher Information decreases per individual 
RSS measurement with increasing sampling density until the 
localization precision stabilizes. When n « n e g, ( |48j ) provides 
practically the same results as the CRLB for independent 
measurements (<2mm difference in all cases). Hence, one can 
assume independent measurements as long as the inverse of 
the sampling rate is large compared to the correlation length. 

The bound on estimator covariance becomes of interest 
when the density of measurements surpasses the bound im¬ 
posed by the spatial coherence region of far-field radiation of 
the transmitters. When the number of measurements goes to 
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infinity, ( |43| reduces to 

lim COV(0) = pa 2 , (49) 

n—f oo 

which does not approach zero. Hence, Bienayme’s theorem 
reveals the link of measurement variance to diffraction, the 
sampling theorem, and Fisher Information. 

G. Bias and Efficiency of Estimator 

We performed 10, 000 simulation runs using the measure¬ 
ment setup described in Fig. [I] to quantify (1) the bias and 
(2) the efficiency of our estimator in ([35}. We performed 
two sets of simulations, one with independent noise and one 
with correlated noise expressed by ([25}. We did not perform 
simulations with cross-correlation function ( [26} , because the 
covariance matrix is singular above the Nyquist sampling 
rate and thus for our measurement setup. Note that these 
simulations do not provide an experimental validation of the 
model used. An estimator is unbiased if © 

E(0) -0 = 0. (50) 

We are interested in the bias of the estimated blind radio 
position (x mle ,y mle ), so we define the bias as 

BIAS(a; mle , t/ mle ) = J (E(a; mle ) - x ) 2 + (E(y mle ) - y ) 2 . 

(51) 

Simulations show that the bias is of the order of 
BIAS(x mle , y mle ) = 0.1mm for independent and correlated 
noise. Estimator efficiency is defined as the difference between 
the covariance of the estimator and CRLB. We quantify this 
difference by 

Estimator efficiency = |i?MS'£ l S i mu i a ted — \/tr(.F - 1 (0))|. 

(52) 

Simulations show that the estimator and CRLB differs at most 
by Estimator efficiency < 1mm for both independent and 
correlated noise. Hence, our estimator in ( [35} has a negligible 
bias and a high efficiency in our simulation environments with 
independent and correlated noise. 

V. EXPERtMENTAL RESULTS 

This section presents the experimental results of the two- 
dimensional measurement setup described in Section [TT] The 
first subsection presents the far-held spatial cross-correlation 
function of the RSS signals and estimate the spread of these 
spatial correlations. We then determine the spatial Fourier 
transform of these cross-correlations to look for an upper 
bound at a spatial frequency of k = 27 t/(Ao/ 2 ), in line with 
the spatial resolution being bounded by the diffraction limit. 
Finally, we determine the global RMSE of ( |40} and show its 
asymptotic behavior to the diffraction limit with the increase 
of the density of sampling points. We compare this experi¬ 
mentally determined RMSE with the RMSE determined by 
the CRLB for independent and correlated noise. We show that 
the CRLB for independent noise underestimates the RMSE 
computed from our measurements. 



Fig. 4. Spatial cross-covariance of RSS measurements as a function of cor¬ 
relation distance in wavelengths. The black curve shows the measured cross¬ 
correlations using all 60 million RSS measurements. The red curve shows 
the measured cross-correlations using 1 RSS measurement per reference radio 
position. 


A. Spatial Correlations and Spatial Frequency Distribution 

Fig-0 shows the spatial cross-covariance function between 
power flows as a function of distance in wavelengths. The spa¬ 
tial cross-covariance function is calculated using the deviations 
from the LNSM expressed by ( |22} and calibrated propagation 
parameters expressed by ( [34} . Hence, spatial cross-covariance 
and thus cross-correlations are distance independent by cross 
correlating the deviations from the calibrated LNSM. Fig. 0 
shows that the spatial cross-correlations go to a minimum over 
a distance of roughly half the mean wavelength, which corre¬ 
sponds to the diffraction limit. The small difference between 
the black and red curve indicates that noise resulting from 
repetitive multiplexed measurements over time is negligible 
compared to the cross-correlated noise in the far field. 

In the case of Fisher Information and thus CRLB analysis, 
information is additive when measurements are independent 
16). One usually assumes signal models with independent 
noise Q7), j2D). Hence, localization precision increases 
with an ever increasing amount of independent measure¬ 
ments over a finite measuring range. However, when measure¬ 
ments are correlated, information and localization precision 
gain decrease with increasing correlations. Fig. 0 shows that 
when space-measurement intervals become as small as the 
diffraction limit, measurements become spatially correlated 
and mutually dependent. Our measurements show that RSS 
signal measurements are spatially correlated over a single¬ 
sided region of roughly half the mean wavelength. It corre¬ 
sponds to the diffraction pattern expressed by © derived 
in Section III-E Therefore, increasing reference radio density 
beyond one per half a wavelength has a negligible influence 
on Fisher information gain and thus on localization precision. 
Our experimental results in the next subsection confirm this. 
In case of independent measurements, Fig. 0 would show an 
infinitely sharp pulse (Dirac delta function). 

Fig. 0 shows the measured power spectrum of cross- 
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Fig. 5. Measured power spectrum of cross-correlated noise in RSS signals 
as given by the spatial Fourier transform of the cross-covariance of Fig. p] 
Theoretically, it is represented by tri(k), as it is the Fourier transform ot 
([26}. The vertical line represents the spatial frequency that corresponds to the 
diffraction limit, which forms an upper bound in the spectrum. 


Fig. 6. Experimentally determined RMSE as a function of all RSS measure¬ 
ments over space (red curve) with the circles indicating RMSE processing 
one RSS measurement per reference radio position. Estimated RMSE using 
Bienayme’s theorem and Fisher Information as in (48) w ith cross-correlations 
equal to the diffraction pattern as expressed by (26) (blue curve), CRLB 
analysis as a function of 1 independent RSS measurement per reference radio 
position (green curve). 


correlated noise in RSS signals, \p(k)\, i.e. the spatial Fourier 
transform of the cross-covariance of Fig. [4] Theoretically, it 
is represented by the spatial Fourier transform of ( |26| . This 
figure shows that the energy is mainly located in lower spatial 
frequencies and it diminishes over a single-sided interval of 
k < 27t/(Ao/2). This upper bound corresponds to the diffrac¬ 
tion limit. The Nyquist sampling rate provides an estimate of 
the minimum sampling rate to fully reconstruct the power-flow 
signal over space without loss of information, which equals the 
single-sided bandwidth of our power spectrum. The vertical 
black line represents the spatial frequency associated with the 
Nyquist sampling rate, which equals k = 27r/(Ao/2). Our 
experimental results in the next subsection confirm this by 
showing that the localization precision does not increase by 
sampling beyond this sampling rate. In case of independent 
measurements. Fig. [5] would show a uniform distribution. 

B. Localization Precision 

Fig. [6] shows the experimental results of our measurement 
setup described in Section [II] Localization precision is given 
as the inverse of RMSE, which is computed from ( |40| . 

The red curve in Fig. [6] shows the RMSE as a function 
of the number of measurements per wavelength. The RMSE 
decreases with increasing number of RSS measurements over 
space until sufficient measurements are available. At a critical 
measurement density, the bound on localization precision be¬ 
comes of interest. The RMSE (red and black curves) converges 
asymptotically to roughly half the mean wavelength as one 
would expect from the diffraction limit. 

The RMSE represented by the red curve is based on process¬ 
ing all 25000 RSS signal measurements per reference radio 
position. The RMSE represented by the black curve is based 
on processing one RSS signal measurement per reference 
radio position. The negligible difference between the red and 
black curves shows that the number of repeated RSS signal 


measurements per reference radio has a negligible influence 
on the RMSE. 

The CRLB for independent measurements starts deviating 
from the RMSE (red and black curves) when the sampling 
rate is increased beyond one RSS signal per half the mean 
wavelength (see black dotted curve), as one would expect 
from the diffraction limit and the Nyquist sampling rate over 
space 114 §2.4], Spatial correlations between RSS signals 


increase rapidly with increasing RSS measurement density 
beyond one sample per half the mean wavelength (Fig. [4|. 
Correlated RSS signals cannot be considered as independent. 
|2J has shown that Fisher information is upper bounded by 
uncertainty in line with ( |2T| >. As o> is upper bounded in spatial 
frequencies by 2/>; () , the spread or uncertainty in position, ay, 
is lower bounded. Hence its inverse is upper bounded, as is 
Fisher Information. Coherence and Speckle theory have shown 
that uncertainty is lower bounded by the diffraction limit. 
Hence, the CRLB at a sampling density of one sample per 
half the mean wavelength (vertical dotted black curve) should 
equal the measured bound on RMSE. We define the measured 
bound on RMSE as the RMSE processing all 60 million 
measurements. Our measurements show that the difference 
between this CRLB and the measured bound on RMSE is 
2-3% (1mm). Hence, our experiments validate the theoretical 
concepts introduced by (2). Our experiments reveal evidence 
that the CRLB cannot be further decreased by increasing the 
number of measurements. On the other hand, at 25 RSS signal 
measurements per wavelength, the RMSE is a factor of four 
higher than the one calculated by the CRLB for independent 
noise. This difference cannot be explained by the difference 
between the covariance of the estimator and the CRLB (see 
Section |IV| >. The difference can, however, be explained by 
calculating the degrees of freedom of our localization setup 
using Bienayme’s theorem and the lower bound on correlation 
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length. When one substitutes this lower bound in the [26) and 
substitutes ([26]) in ( [38] ), one obtains the asymptotic behavior 
of the CRLB with correlated noise as is shown in Fig. [3] The 
measured asymptote in Fig. |6]deviates 2-3% from Bienayme’s 
theorem and from the CRLB for correlated noise. 

All 60 million power measurements and Matlab files are 
arranged in a database at Linkoping University m- 


VI. Discussion 


Our novel localization setup of using reference radios on the 
circumference of a two-dimensional localization area instead 
of setting them up in a two-dimensional array worked well. 
This implies that one does not have to know the phases and 
amplitudes on the closed surfaces around extended sources 
to reconstruct the positions of these extended sources. Our 
experiments show that it suffices to measure time-averaged 
power flows when a localization precision of about half a 
wavelength is required (in our case 6 cm). Time averaging on 
a time scale that is long compared to the temporal coherence 
time is usually employed in RSS localization, so that phase 
information is lost. We expect such a setup to work well when 
all radio positions are in LOS. For a practical implementation 
in NLOS environments, we refer to 123 §5]. The propagation 


equations of Section [ill] are than applied to design an efficient 
setup and algorithm to locate the blind radios. 


In Section III we distinguished between primary and sec¬ 
ondary cross-correlated noise in the far field. In our ex¬ 
perimental setup, where signal levels are large compared to 
noise levels, it is reasonable to neglect thermal and quantum 
noise. Scattering from the surface roughness of the primary 
source cannot provide the the spatial frequency bandwidth of 
Fig. [5] because of the small size of the transmitting antenna. 
We expect that the cross-correlated noise originates from the 
surface roughness of the large area between the transmitter and 
receiver. Hence, our 2400 time-averaged RSS measurements 
are correlated over space in line with the derived correlation 
model of ( [26] ) as was verified by our experiments. 

Noise can originate from a variety of sources. In our 
setup, the noise level of the antenna plus electronics on a 
typical 802.15.4 radio is about 35dBm below signal level as 
specified in ]30[ . The reproducible ripples on our RSS signals 
originated in part from small interference effects from undue 
reflections from mostly hidden metallic obstacles in the mea¬ 
suring chamber. Such obstacles like reinforced concrete pillars 
could not be removed. In a real indoor office environment, 
these interference effects usually dominate RSS signals. In 
these environments the cross-correlated noise in RSS signals 
is still determined by the diffraction limit, which in turn, 
determines the bound on localization precision. Multi-path 
effects and fading do not change the stochastic properties of 
the secondary extended sources where the noise is generated. 


As explained in the literature |18|, such models are usually 


clarified by ray-tracing. Ray-tracing implies the geometrical- 
optics approximation, so that the wave lenghth is set to zero 
and diffraction effects are not considered. When we set the 
correlation length of the exponentially decaying correlations 
as expressed by |4j equal to the diffraction limit, the CRLB 


converges to practically the same bound on localization preci¬ 
sion. It is remarkable that this bound is revealed by our simple 
experiment using various noise and signal models. 

Our experiments reveal evidence that the diffraction limit 
determines the (1) bound on localization precision, (2) sam¬ 
pling rate that provides optimal signal reconstruction, (3) size 
of the coherence region and uncertainty, (4) upper bound on 
spatial frequencies and (5) upper bound on Fisher Information 
and lower bound on the CRLB. With the mathematical work 
of |2J, a rigorous link between Fisher Information, CRLB 
and Uncertainty was established. As uncertainty has a lower 
bound according to coherence and speckle theory, so must 
Fisher Information have an upper bound, so that the CRLB 
has a lower bound, all when the noise processes are assumed 
to be Gaussian distributed. Our experiments were able to 
validate those theoretical concepts. This is further confirmed 
by applying Bienayme’s theorem. An interesting observation 
is that despite the fact that power flows have twice the spatial 
bandwidth than radiation flows determined by amplitudes, 
their spatial coherence regions are the same in size, as are 
their degrees of freedom. 

Our paper describes a novel experimental and theoretical 
framework for estimating the lower bound on uncertainty for 
localization setups based on classical wave theory without 
any other prior knowledge. Bienayme’s theorem and existing 
CRLB for signals with correlated Gaussian noise reveal our 
postulate [48) , so that the lower bound on uncertainty corre¬ 
sponds to the upper bound on Fisher Information. Our exper¬ 
imental results cannot be explained by existing propagation 
models with independent noise. It took almost three weeks 
in throughput time to perform the 60 million measurements, 
generating 2400 x 50 x 500 repetitive multiplexed measure¬ 
ments. We tried to minimize multi-path effects by avoiding the 
interfering influence of ground and ceiling reflections. Making 
sure to minimize reflections of other metallic obstacles in the 
measuring chamber were challenging but could be overcome. 
This allowed us to reveal a performance bound in theory and 
experiment for measurements without phase information. 

VII. Conclusion 

Our novel two-dimensional localization setup where we 
positioned the reference radios on the circumference of the 
localization area rather than spreading them out in a two- 
dimensional array over this area worked well in our LOS 
setup. Our measurements show that localization performance 
does not increase indefinitely with an ever increasing number 
of multiplexed RSS signals (space and time), but is limited 
by the spatial correlations of the far-field RSS signals. When 
sufficient measurements are available to minimize the influ¬ 
ence of measurement noise on localization performance, the 
bound on localization precision is revealed as the region of 
spatially correlated far-field radiation noise. The determination 
of this region of spatial correlations is straightforward and 
can be directly calculated from RSS signals. Within this 
region of correlated RSS signals, assumptions of independent 
measurements are invalidated, so that the CRLB for indepen¬ 
dent noise underestimates the bounds on radio localization 
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precision. This underestimation is removed by accounting for 
the correlations in the noise. The bound on the correlations 
is given by the fundamental bound on the correlation length 
that we derived from first principles from the propagation 
equations. The CRLB is linked to the uncertainty principle 
as measurement variance is directly related to this principle as 
we showed. The bounds on the precision of RSS- and TOF- 
based localization are expected to be equal and of the order 
of half a wavelength of the radiation as can be concluded 
from our experiments and underlying theoretical modeling. 
Sampling beyond the diffraction limit or the Nyquist sampling 
rate does not further resolve the oscillation period unless near- 
instantaneous measurements are performed with the a priori 
knowledge that signal processing is assumed to be based on 
non-linear mixing. Future research is aimed at the inclusion 
of strong interference effects such as show up in practically 
any indoor environment. 
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